Computational Investigations of the Primary Excited States of Poly(para-phenylene 

vinylene) 



in 
o 
o 

(N 

o 

Q 

(N 



D 
i 

CO 

c3 



I 

a 

o 

o 



> 

On 
(N 

in 
o 



a 

o 

o 



X 



Robert J. Bursill 1 and William Barford 2 
1 School of Physics, University of New South Wales, Sydney, New South Wales 2052, Australia 
1 Department of Physics and Astronomy, University of Sheffield, Sheffield, S3 7RH, United Kingdom 

The Pariser-Parr-Pople model of 7r-conjugated electrons is solved by the density matrix renormal- 
ization group method for the light emitting polymer, poly(para-phenylene vinylene). The energies of 
the primary excited states are calculated. When solid state screening is incorporated into the model 
parameters there is excellent agreement between theory and experiment, enabling an identification 
of the origin of the key spectroscopic features. 

PACS numbers: 71.10.-w, 71.20.Rv, 71.35.-y, 78.67.-n 



I. INTRODUCTION 

Since the discovery of electroluminescence in 
poly(para-phenylene vinylene) (PPV) in 1990i, a 
variety of experimental and theoretical techniques have 
been deployed to investigate the physics of the primary 
photoexcitations of the phenyl-based light emitting 
polymers. Linear and non-linear optical spectroscopies 
have revealed the energies and symmetries of the 
dominant dipole-allowed transitions, as well as the 
important dipole-forbidden transitions that participate 
in non-linear optical processes. 

All phenyl-based light emitting polymers exhibit a 
characteristic absorption spectra. These show two dom- 
inant excitations polarized parallel to the chain axis (at 
2.8 eV and 6.1 eV in PPV-). There are also two in- 
termediate weaker transitions (at 3.6 eV and 4.8 eV in 
PPV 2 ). The higher of these intermediate transitions is 
predominantly polarized perpendicular to the chain axis, 
whereas the polarization of the lower transition is less 
well-defined. (In PPV it is predominately polarized par- 
allel to the long-axis^.) Furthermore, the strength of this 
lower transition is enhanced by chemical substitution. 
Electroabsorption2> and two-photon absorption— reveal a 
dipole-forbidden state at approximately 0.7 eV above the 
lowest dipole-allowed transition. These two states are 
usually labelled the rn^Ag and \ ^B U states, respectively. 
Finally, a triplet state has been observed at 0.7 eV be- 
low the \ X B U excitation (see reS for references), with 
another dipole connected triplet state 1.4 eV higher in 
energy 6 . (Notice that this higher triplet state, labelled 
as the m 3 A g state, is virtually degenerate with its singlet 
counterpart, namely the m}A g state. These two states 
are often referred to as the singlet and triplet charge- 
transfer states. Their degeneracy can be explained by the 
fact that they are the lowest pseudomomentum branches 
of the n = 2 Mott-Wannier excitons, which have odd par- 
ity electron-hole wavefunctions and therefore experience 
no exchange interactional.) 

An important early insight into the nature of the pri- 
mary photoexcitations of the phenyl-based systems was 
provided by Rice and Gartstein&ii, who argued that they 
can essentially be understood as arising from the dereal- 



ization of the primitive benzene excitations. Kirova and 
Brazovskii, on the other hand, have argued that a conven- 
tional semiconductor band picture of bound particle-hole 
excitations is a more appropriate description 10 . Other 
theoretical work on PPV includes, a single configuration 
interaction (CI) calculation of the Pariser-Parr-Pople 
model by Chandross and Mazumdar 11 ; density matrix 
renormalization group calculations on reduced molecular- 
orbital models by Barford et alH; quantum chemistry 
calculations on the INDO Hamiltonian by Beljonne et 
ali&, and Weibel and Yaroni4; and ah initio Bethe- 
Salpctcr equation (BSE) calculations by Rohlfing and 
Louiei^. These theoretical predictions will be discussed 
more fully in Section IlIII when we discuss the results of 
the calculations presented here and their comparisons to 
experiment. 

In this paper we present density matrix renormaliza- 
tion group (DMRG) calculations on the full Pariser-Parr- 
Pople model of PPV. The DMRG method for calculating 
the full electronic spectra of conjugated polymers has a 
number of advantages over its competitors. It is more 
accurate than single or double CI calculations, and does 
not suffer from finite-size consistency errors; it makes 
no assumptions about the nature of the excitations (un- 
like the BSE method, which assumes that the excita- 
tions are particle-hole pairs), thus it is able to accurately 
model highly correlated excited states; and unlike stan- 
dard quantum chemistry techniques, it is able to calcu- 
late the excited spectra of large molecules. The DMRG 
method is most readily suited to reduced basis Hamilto- 
nians, such as the 7r-electron Pariser-Parr-Pople model. 
This is not necessarily a disadvantage over ab initio meth- 
ods, as when correctly parameterized the Pariser-Parr- 
Pople model makes very accurate predictionsii*^. More- 
over, the Pariser-Parr-Pople model possess particle-hole 
symmetry, which means that spatial, particle-hole and 
spin-flip symmetries can be employed to target high-lying 
excited states in the DMRG method (see refiii, for exam- 
pie). 

DMRG calculations on the full Pariser-Parr-Pople 
model were presented for poly(para-phenylene) by Bur- 
sill and Barford 18 . This paper extends that approach 
to PPV. In the next section we define the Pariser-Parr- 
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FIG. 1: The band structure of PPV with all nearest neighbor 
bond integrals, t, equal. The low-lying particle-hole transi- 
tions are labelled 1, 2 and 3. 

Poplc model and briefly describe our implementation of 
the DMRG method. In Section ITTT1 we describe and dis- 
cuss our results, and compare them to other approaches, 
concluding in section llVl 

II. MODEL AND METHODOLOGY 

A. The Pariser-Parr-Pople Model 

The Pariser-Parr-Pople model is a 7r-electron model of 
conjugated polymers, defined by, 

H = — tij [cjjj.CjCT + CjQ.Ciu 

<ij> cr 

+ U J2(n n -l/2) (nn-1/2) 

i 

+ lj2 v iM- 1)^-1), (i) 

where <> represents nearest neighbors, Ci„ destroys a 
7r-electron on site i, rii a — c\ a Ci„, and rii = n.^ + Ujj,. 

We use the Ohno parameterization for the Coulomb 
interaction, defined by, 

Va = C// v /l + (C/er. iJ /14.397) 2 , (2) 

where is the inter-atomic distance (in A) , U is the on- 
site Coulomb interaction (in eV), and e is the dielectric 
constant. This interaction is an interpolation between an 
on-site Coulomb repulsion, U , and a Coulomb potential, 
e 2 /4ireeorij as — » oo. 

The band structure, obtained in the non-interacting 
limit (U — 0), is shown in Fig. ^ The pair of non- 
bonding bands are a consequence of the Z?2?i symmetry 



of the Hamiltonian in the non-interacting limit with only 
nearest- neighbor bond integrals. The low- lying particle- 
hole excitations are labelled 1 (for transitions involving 
di and d* ) , the degenerate pair labelled 2 (for transitions 
involving d\ and I*, and I and d\) 1 and 3 (for transi- 
tions involving I and I*). Coulomb interactions lift the 
degeneracy of the intermediate pair (which become the 
intermediate pair of transitions described in Section [IJ , 
while the transitions 1 and 3 become the low and high 
energy transitions described in Section [I] Evidently, to 
obtain a more realistic description of the excited states 
it is necessary to include Coulomb interactions. 

In the following we use two sets of parameters to 
model the excited states. The first set, called the op- 
timized parameters, were derived by fitting the predicted 
Pariser-Parr-Pople excitation energies of stilbene to the 
experimental spectrum of stilbene in vacuo<i2i. These are, 
t p = 2.539 eV, t d = 2.684 eV, t s = 2.22 eV, U = 10.06 
eV, and e = 1 (where the phenyl, double and single bond 
integrals are defined in Fig. [2J . Although giving instruc- 
tive results, this parameter set fails to quantitatively pre- 
dict the excitation energies of polymers in the solid state, 
as they fail to account for the solvation effects of the sur- 
rounding dielectrio22i2i. The other parameter set, called 
the screened parameters, were derived by Chandross and 
Mazumdar— to account for solid state solvation effects. 
These parameters are, t p = 2.4 eV, td — 2.6 eV, t s = 2.2 
eV, U = 8 eV, and e = 2. 



B. The Density Matrix Renormalization Group 
(DMRG) Method 

Eq. 0} is solved by the DMRG method for chains of up 
to 28 phenyl rings (i.e. 222 sites). The DMRG method 
is an efficient truncation procedure for solving quantum 
lattice Hamiltonians, especially in one-dimension^. The 
details of the DMRG implementation for this problem, 
particularly the procedure of deriving an optimized re- 
duced basis for the phenyl rin g 18 - 23 - 24 ; are described in 
detail in re£i&. Also shown in re£i& are comprehensive 
DMRG convergence tests that are applicable to phcnyl- 
based systems. 

As described in Section[IJ the Pariser-Parr-Pople model 
possess spin-flip and particle-hole symmetries, which we 
exploit in the DMRG method to target excited states. 
With only onsite and nearest neighbor Coulomb interac- 
tions the Pariser-Parr-Pople model applied to the PPV 
structure formally possess T)ik symmetry. In this limit 
this means that the eigenstates are labelled with the 
spatial symmetry assignments, A g , B\ Ul Biu, and £>3 g . 
Longer range Coulomb interactions with the PPV struc- 
ture reduce the Dih symmetry to Ci symmetry, and thus 
the true spatial symmetry assignments are A g and B u . 
However, it is computationally expedient, with very little 
loss of accuracy, to retain D2/1 symmetry while incorpo- 
rating long range Coulomb interactions. This is achieved 
by taking the bond angle between the single and double 
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FIG. 2: The DMRG calculated transition energies of para- 
phenylene vinylene oligomers as a function of inverse chain 
length, calculated from the Pariser-Parr-Pople model with op- 
timized parameters: U = 10.06 eV, t p = 2.539 eV, U = 2.684 
eV, t s — 2.22 eV, and dielectric constant, e = 1. The excited 
states are identified with the spectroscopic features shown in 
Fig. 3 of re£ and Fig. 1 of rei». The symmetry assignments 
of the states shown in brackets are the symmetries the eigen- 
states would have if PPV had D^h rather than Ci symmetry. 
The inset shows the oligo-phenylene vinylene structure with 
the bond integrals t p , td and t s . 



bonds in the vinylene unit to be 180°, rather than 120°. 
To ensure that the overall molecular size is consistent 
to the original structure the single and double lengths 
are reduced to 1.283 A and 1.194 A, respectively. The 
phenyl bond length is 1.40 A. The next section describes 
the calculated results. 



III. RESULTS AND DISCUSSIONS 
A. DMRG Calculations 

Fig. [5] shows the DMRG calculated excitation ener- 
gies of oligo(para-phcnylcnc vinylcncs) using the Pariser- 
Parr-Pople model with unscreened parameters {U — 
10.06 eV, t p = 2.539 eV, t d = 2.684 eV, t s = 2.22 eV, 
and e = 1). (N = 2 corresponds to stilbene.) The 1 X B~ 
state is the strong lowest energy dipole allowed transition 
(labelled I in Fig. 3 of refi). The strong transition dipole 
moment between this state and the 2 1 J 4+ state indicates 
that the 2 1 ^4+ is the state labelled 'mlAg' in non-linear 



optical spectroscopies (see, for example, Fig. 1(a) in refi). 
The energies of the 1 X B~ and 2 1 ^4+ states initially reduce 
rapidly as a function of chain length, indicating that these 
states reduce their kinetic energy by readily delocalizing 
the exciton wavefunction along the chain. This is pri- 
marily because these states are largely constructed from 
particle-hole transitions between the bonding and anti- 
bonding bands, labelled c?i and d\ in Fig.^ 5 . A descrip- 
tion of these states as strongly bound band exciton states 
is therefore appropriate in the long-chain limit. Indeed, 
the ^ l B~ state is the lowest pseudomomentum branch 
of the n = 1 Mott-Wannier exciton, while the 2 : A+ is 
the lowest pseudomomentum branch of the n — 2 Mott- 
Wannier excitonLii. (Higher lying pseudo-momentum 
branches of the n = 1 and n = 2 excitons for DMRG cal- 
culations on poly (par a-phenylene) are illusrated in Fig. 
5 of re£i^.) Notice that this assignment places a lower 
bound on the binding energy of the n = 1 exciton as 
E(mA g ) - E(1B U ). 

We next consider the state labelled II and shown by 
diamonds in Fig. [21 This state has 1 B+ symmetry, but 
would have X B^ U symmetry if PPV had D-xh symmetry. 
From a band picture analysis, this state is composed of 
an antisymmetric combination of particle- hole transitions 
from e?i to I* and I tod*. At the Pariser-Parr-Pople model 
level, it has positive particle-hole symmetry, and is there- 
fore dipole-forbidden. However, it is weakly dipole al- 
lowed in real systems, and its oscillator strength is further 
enhanced by chemical substitution. Its excitation energy 
is virtually independent of chain length, because, first its 
wavefunction is composed of non-bonding orbitals, and 
therefore there is no hybridization via one-electron trans- 
fer terms, and second since it has a very small oscillator 
strength, resonant exciton transfer along the chain is vir- 
tually inoperative. These two features (weak oscillator 
strength and a chain independent energy) strongly indi- 
cate that the second spectroscopic feature (labelled II in 
Fig. 3 of refi) originates from this state. 

The state labelled III and shown by triangles in Fig. 
El has 1 B~ symmetry, but would have 1 B^ U symmetry if 
PPV had D^h symmetry. From a band picture analysis, 
this state is composed of a symmetric combination of 
particle-hole transitions from d\ to I* and I to d\. It 
has negative particle-hole symmetry, and therefore has 
an allowed dipole transition from the ground state. This 
state is the third spectroscopic feature of PPV (labelled 
III in Fig. 3 of re£). 

The state labelled IV and shown by crosses in Fig. [21 
has l B~ symmetry, but would have 1 B^ A1 symmetry if 
PPV had D2h symmetry. From a band picture analysis, 
this state is composed of particle-hole transitions from 
I to I*. It has negative particle- hole symmetry, has a 
strong dipole transition from the ground state, and is 
polarized along the chain axis. This state is the fourth 
spectroscopic feature of PPV (labelled IV in Fig. 3 of 
refi), as is usually referred to as the 'Frenkel' exciton, as 
its particle-hole wavefunction is typically confined to a 
single phenyl-ring. 
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FIG. 3: The DMRG calculated transition energies of para- 
phenylene vinylene oligomers as a function of inverse chain 
length, calculated from the Pariser-Parr-Pople model with 
screened parameters: U = 8 eV, t p = 2.4 eV, t d — 2.6 eV, 
i s = 2.2 eV, and e = 2. The symbols are defined in the inset 
of Fig. |2l The excited states are identified with the spectro- 
scopic features shown in Fig. 3 of re£ and Fig. 1 of reA 



Finally, the two triplet states, 1 3 B+ and m 3 A~, are 
also shown in Fig. [21 

The results shown in Fig. |3 are obtained from the 
Pariser-Parr-Pople model with optimized parameters 
that are not parametrized to model solid state solvation 
effects. The dielectric response of the environment is pre- 
dicted to red-shift the intra-molecular excitations of a 
single chain by up to 0.1 eV for the l 1 ^ state, 0.7 eV 
for the m}A g state, and 1.5 eV for the charge gap2£L2I. 
We would therefore not expect these calculations on a 
single chain to compare directly with experimental ob- 
servations. 

Chandross and Mazumdar attempted to model solva- 
tion affects by adjusting the Pariser-Parr-Pople param- 
eters and introducing a static dielectric constantsAi. Al- 
though the introduction of a static dielectric constant is 
not a faithful representation of the dielectric response 
(because the timescale for the particle-hole motion is 
very similar to the timescale of the dielectric response^) , 
these screened parameters do lead to Pariser-Parr-Pople 
model predictions that are remarkably close to the ex- 
perimental observations. 

Fig. |21 shows the DMRG calculated results with the 
screened parameters ({7 = 8 eV, t p = 2.4 eV, f<j = 2.6 eV, 
t s = 2.2 eV, and e = 2). Comparing this figure with Fig. 
|21we see that the ^B~ and 2 1 A+ states are red-shifted 
by ca. 0.1 eV and 0.6 eV, respectively, as expected from 
solvation effects. The excited states are identified with 
the spectroscopic features shown in Fig. 3 of re£. Their 
energies are remarkably consistent with the experimental 
values of 2.8, 3.6, 4.8, and 6.1 eV obtained in reS, lend- 
ing additional credence to our assignments of the excited 



state origins of the features, discussed above. 



B. Other Approaches 

We now compare our DMRG predictions of the Pariser- 
Par-Pople model with other approaches. Our results us- 
ing the screened parameters are consistent with those 
of Chandross and Mazumdar—. They used SCI on the 
Pariser-Par-Pople model. For an eight-unit oligomer they 
calculate the l 1 ^ state at 2.7 eV, an w}A g state at 3.3 
eV, and the n}B u state at 3.6 eV. The 1 X B U and m 1 A g 
states are the n = 1 and n — 2 excitons, while the n l B u 
state coincides with the charge-gap and therefore indi- 
cates the onset of the particle-hole continuum. They also 
predict the l 3 B u state at 1.4 eV. 

An ab initio BSE calculation by Rohlfing and LouieA^ 
on a PPV polymer predicts dipole allowed and forbid- 
den singlet excitons at 2.4 eV and 2.8 eV, respectively, 
with the quasi-particle gap at 3.3 eV. They also predict 
triplet excitons at 1.5 eV and 2.7 eV. The 2.4 eV and 
2.8 eV singlet excitons are the ~\}B U and m 1 A g states, 
respectively, while the 1.5 eV and 2.7 eV triplet excitons 
are the 1 3 B U and m 3 A g states, respectively. The m 1 A g 
and m 3 A g states are nearly degenerate, as predicted by 
the Mott-Wannier exciton theory for odd parity particle- 
hole wavefunctions. Using the same technique with a 
screened electron-hole interaction van der Horst et aZ— 
predict a l^Bu binding energy in PPV of 0.48 eV. 

The origin of the higher-lying transitions has also 
been investigated. Rohfiing and Louiei^, and Weibel 
and Yaroni4 predict that peak II in PPV arises from 
an exciton caused predominately by the antisymmet- 
ric combination of particle-hole transitions from d\ to 
I* and I to d\. Weibel and Yaroni4 have also investi- 
gated the effects of breaking particle-hole symmetry on 
the oscillator strength and polarization of peak II. Using 
the semiempirical INDO Hamiltonian on nonplanar di- 
hydroxy-PPV, their calculations indicate that chemical 
substitution and mixing of the 7r and a orbitals enhances 
the oscillator strength, as originally suggested by Gart- 
stein et al& Moreover, as illustrated in Fig. 5 of refii^, this 
peak becomes predominately polarized along the chain 
axis, in agreement with experiment. 



IV. CONCLUSIONS 

Gathering together the various theoretical predictions 
of the origins of the primary excited states of PPV we 
now interpret the key spectroscopic features of PPV. 

Peak I corresponds to the low-energy dipole active 
1 X B~ state. This is the lowest pseudomomentum branch 
of the family of n = 1 Mott-Wannier singlet exci- 
tons resulting from the Coulomb attraction between the 
particle- hole excitation from the valence (d±) to the con- 
duction (d*) bands. Approximately 0.7 eV higher in en- 
ergy is the m^^Ag state, identified by electroabsorption 2 , 
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two-photon absorption and photoinduced absorption 4 . 
The Pariser-Parr-Pople model calculations described in 
this paper suggest that this state is the 2M+ state, which 
is the lowest pseudomomcntum branch of the family of 
n = 2 Mott-Wannier excitons. This assignment places a 
lower bound on the spectroscopically determined binding 
energy of the 1 X B U exciton of 0.7 eV. Approximately 0.7 
eV below the 1 X B~ exciton is the 1 3 B+ triplet, indicat- 
ing a large exchange energy characteristic of correlated 
states. This state is the lowest pseudomomentum branch 
of the family of n = 1 Mott-Wannier triplet excitons. 
Photo-induced absorption from the l 3 B+ triplet indi- 
cates another triplet, the l 3 ^4g state, at approximately 
1.4 eV higher in energy, and essentially degenerate with 
the 2 1 A+ state. This triplet state is the lowest pseudo- 
momentum branch of the family of n — 2 Mott-Wannier 
triplet excitons. As expected from Mott-Wannier exciton 
theory in one-dimensiorAii, the odd particle-hole parity 
singlet and triplet (charge-transfer) excitons are virtually 
degenerate. The n}B u state at 0.1 eV higher in energy 
than the m}A g state in PPV— indicates binding energies 
of ~ 0.8 eV and 0.1 eV for the n = 1 and n = 2 singlet 
excitons, respectively. Higher in energy are the excita- 



tions associated with peaks II and III. These arise from 
antisymmetric and symmetric combinations of particle- 
hole transitions from d\ to I* and I to d*, respectively. 
Finally, peak IV is the intraphenyl, or Frenkel, exciton. 
This identification of the primary excited states of PPV 
also applies to other phenyl-based systems, e.g. PPRiLi&. 

The calculations presented in this paper have entirely 
concerned vertical transitions, where the excited states 
are calculated in the geometry of the ground state. Of 
course, electron-lattice relaxation is an important process 
in polymers that determines Stokes shifts, self-trapping, 
and exciton and charge transfer between polymers. The 
consequences of electron-lattice relaxation on the excited 
states energies and structures have been considered in 
refs2£. 
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